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Detailed analyses of the WMAP data indicate possible oscillatory features in the primordial 
curvature perturbation, which moreover appears to be suppressed beyond the present Hubble radius. 
Such deviations from the usual inflationary expectation of an approximately Harrison-Zeldovich 
spectrum are expected in the supergravity-based 'multiple inflation' model wherein phase transitions 
during inflation induce sudden changes in the mass of the inflaton, thus interrupting its slow-roll. In 
a previous paper we calculated the resulting curvature perturbation and showed how the oscillations 
. arise. Here we perform a Markov Chain Monte Carlo fitting exercise using the 3-year WMAP data 

' to determine how the fitted cosmological parameters vary when such a primordial spectrum is used 

' as an input, rather than the usually assumed power- law spectrum. The 'concordance' ACDM model 

is still a good fit when there is just a 'step' in the spectrum. However if there is a 'bump' in the 
^ ' spectrum (due e.g. to two phase transitions in rapid succession), the precision CMB data can be 

Q ' well-fitted by a flat Einstein-de Sitter cosmology without dark energy. This however requires the 

, Hubble constant to be ~ 0.44 which is lower than the locally measured value. To fit the SDSS 

■ data on the power spectrum of galaxy clustering requires a ~ 10% component of hot dark matter, as 

fT^ ' would naturally be provided by 3 species of neutrinos of mass ~ 0.5 eV. This CHDM model cannot 

however fit the position of the baryon acoustic peak in the LRG redshift two-point correlation 
function. It may be possible to overcome these difficulties in an inhomogeneous Lemaitre-Tolman- 
Bondi cosmological model with a local void, which can potentially also account for the SN la Hubble 
diagram without invoking cosmic acceleration. 



o 

C/3 , I. INTRODUCTION 

Precision measurements of CMB anisotropies by the Wilkinson Microwave Anisotropy Probe ( WMAP) are widely 
' accepted to have firmly established the 'concordance' ACDM model — a flat universe with Oa ~ 0.7, fini — 0.3, 
^ . /i ~ 0.7, seeded by a nearly scale-invariant power-law spectrum of adiabatic density fluctuations However the 

model fit to the data is surprisingly poor. For the WMAP-1 TT spectrum, Xoff/^ — 974/893 so formally the 
. ACDM model is ruled out at 97% c.l. Visually the most striking discrepancies are at low multipoles where the lack of 
(SJ ' power in the Sachs- Wolfe plateau and the absence of the expected A-induced late integrated-Sachs- Wolfe effect have 
' drawn much attention. However because cosmic variance and uncertainties in the foreground subtraction are high 

■ on such scales, it has been argued that the observed low quadrupole (and octupole) are not particularly unlikely, see 
J — ' e.g. refs.[l,[3, The excess in fact originates mainly from sharp features or 'glitches' in the power spectrum that 

, the model is unable to fit [H, Although these glitches are less pronounced in the 3-year data release, they are still 
present p]. Hence although the fit to the concordance ACDM model has improved with Xcs/'^ — 1049/982 for the 
WMAP-3 TT spectrum this model still has only a 6.8% chance of being a correct description of the data. This 

■ is less than reassuring given the significance of such a tiny cosmological constant (more generally, 'dark energy') for 
J-j \ both cosmology and fundamental physics. 

The WMAP team state: "In the absence of an established theoretical framework in which to interpret these glitches 
(beyond the Gaussian, random phase paradigm), they will likely remain curiosities" @. However it had been noted 
earlier by the WMAP team themselves Q that these may correspond to sharp features in the spectrum of the 
underlying primordial curvature perturbation, arising due to sudden changes in the mass of the inflaton in the 
'multiple inflation' model proposed in ref.[^. This generates characteristic localized oscillations in the spectrum, as 
was demonstrated numerically in a toy model of a 'chaotic' inflationary potential having a 'step' starting at ^stop with 

amplitude and gradient determined by the parameters Campi and dgiad^ V{(t>) — ^m'^cf)'^ I + Campi tanh ^ ^ ^ 

0. It was found that the fit to the WMAP-1 data improves significantly (by Ax^ = 10) for the model parameters 
^,t3p = 15.5 Afp, Campi = 9.1 X 10-* and dg,ad = 1.4 x lO'^Mp, where Mp = (87rGN)"^/^ ^ 2.44 x 10^^ GeV 0. 
However the cosmological model parameters were held fixed (at their concordance model values) in this exercise. 
Recently this analysis has been revisited using the WMAP-3 data flQ|; these authors also consider departures from 
the concordance model and conclude that there are virtually no degeneracies of cosmological parameters with the 
modelling of the spectral feature [llj . 

However m in the toy model above is not the mass of the inflaton — in fact in all such chaotic inflation models 
with V (x (j)" where inflation occurs at field values 0infl > Mp, the leading term in a Taylor expansion of the potential 
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around (piaR is always linear since this is not a point of symmetry [12| |. The effect of a change in the inflaton mass 
in multiple inflation can be sensibly modelled only in a 'new' inflation model where inflation occurs at field values 
(j> << Mp\ and an effective field theory description of the inflaton potential is possible. The 'slow-roll' conditions are 
violated when the inflaton mass changes suddenly due to its (gravitational) coupling to 'flat direction' fields which 
undergo thermal phase transitions as the universe cools during inflation ,8]. The resulting effect on the spectrum of 
the curvature perturbation was found by analytic solution of the governing equations to correspond to an upward 
step followed by rapidly damped oscillations [13:].^ 

One can ask if spectral features are seen when one attempts to recover the primordial perturbation spectrum directly 
from the data. Such attempts may be divided into two classes. Usually the curvature perturbation, V-R,{k), is given 
a simple parameterisation and fitted to the data, together with the background cosmology, using MCMC likelihood 
analysis. The spectrum has been described using bins in wave number k [l^, [13, [U [22|, wavelets [1^, smoothing 
splines [13] and principal components [1^. However for MCMC analysis to be feasible, the number of parameters 
must be limited, so the reconstructed spectrum has too low a resolution to reveal anything interesting. By contrast, 
'non-parametric' methods assume the background cosmology (usually the concordance ACDM model) so that the 
transfer function is known, and then invert the data to find the primordial perturbation spectrum. Methods that 
have been used are the Richardson-Lucy deconvolution algorithm jJl, H^l, an iterative, semi-analytic process [28Ll29j. 
and a smoothed, least-squares procedure [s^, [3l| . The perturbation spectra produced by the first and third methods 
have a prominent step followed by bumps which are rather reminiscent of decaying oscillations (see Fig. 4 in ref. [l3|). 
These features correspond in fact to the depressed quadrupole and the glitches at low multipole £. 

Encouraged by this we consider possible variations of the primordial perturbation spectrum beyond the limited 
set considered so far, and motivated by the multiple inflation model. We have shown earlier that a phase transition 
in a 'flat direction' field during infiation (in a supergravity framework) generates a step in the primordial spectrum 
followed by damped oscillations [13| . In a physical supergravity model there are many such fiat directions and these 
will undergo phase transitions in rapid succession after the first ~ 10 — 15 e-folds of inflation (assuming they all start 
from the origin due e.g. to thermal initial conditions) Hence we also consider a possible 'bump' in the spectrum due 
to two phase transitions in rapid succession, which raise and then lower the inflaton mass. We wish to emphasise that 
there may well be other physical frameworks wherein one can expect similar features in the primordial spectrum. Our 
intention here is to use a definite and calculable framework, in order to illustrate how the extraction of cosmological 
parameters is dependent on the assumed form of the primordial spectrum. We find that the precision WMAP data 
can be fitted just as well without invoking dark energy if there is indeed a bump in the primordial spectrum around 
the position of the first acoustic peak in the angular power spectrum of the CMB. To fit the position of the first peak 
(assuming a fiat model as motivated by inflation) requires however a low Hubble constant, h ~ 0.44, in contrast to the 
value of h = G.72±0.08 measured by the Hubble Key Project (HKP) in our local neighbourhood [s^l- Although a pure 
CDM model is well known to suffer from excess power on small scales, data from the Sloane Digital Sky Survey {SDSS) 
on galaxy clustering [33| can also be well fitted if there is a 10% component of hot dark matter, as has been noted 
already using the 2dFGRS data [ll, [11]. This is encouraging given the evidence for neutrino mass from oscillations, 
and the required value of ^ 0.5 eV per species is well within the present experimental upper bound of 2.3 eV [36| . 
Such a cold -I- hot dark matter (CHDM) model with a low Hubble constant passes all the usual cosmological tests 
(e.g. cluster baryon fraction and as from clusters and weak lensing [s^) but has difficulty [sj matching the position 
of the 'baryon acoustic oscillation' (BAO) peak observed in the redshift two-point correlation function of luminous 
red galaxies (LRG) in SDSS J38]. We confirm that this is indeed the case but wish to draw attention to the possibility 
that this difficulty may be solved in an inhomogeneous Lemaitre-Tolman-Bondi (LTB) cosmological model wherein 
we are located in an underdense void which is expanding faster than the global rate 13911 . Such a model may also 
account for the SN la Hubble diagram without invoking cosmic acceleration [13, [HI, |42|. |43|. [3, [4^ [4^ . 

II. MULTIPLE INFLATION 

In previous work we have discussed the effective potential during inflation driven by a scalar field in = 1 
supergravity, which has couplings to other flat direction fields having gauge and/or Yukawa couplings These fields 
acquire a potential due to supersymmetry breaking by the large vacuum energy driving infiation and evolve rapidly 



A 'hybrid' infiation model wherein the inflaton is coupled to a 'curvaton' field also yields oscillations together with suppressed power 
on large scales [14j] . A similar phenomenon had been noted earlier for the case where the inflaton potential has a jump in its slope lisll : 
however such a discontinuity has no physical interpretation. The WMAP glitches have also been interpreted as due to the effects of 
'trans-Planckian' physics Il6l ITtII and due to resonant particle production ll^ . 
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to their minima, which are fixed by the non-renormaUsable terms which Hft their potential at large field values. The 
infiaton's own mass thus jumps as these fields suddenly acquire large vacuum expectation values (vevs), after having 
been trapped at the origin through their coupling to the thermal background for the first ~ 10—15 e-folds of inflation. 
Damped oscillations are also induced in the infiaton mass as the coupled fields oscillate in their minima losing energy 
mainly due to the rapid inflationary expansion. The resulting curvature perturbation was calculated in our previous 
work [13] and will be used in this paper as an input for cosmological parameter extraction using the WMAP-3 data. 

In order to produce observable effects in the CMB or large-scale structure the phase transition(s) must take place as 
cosmologically relevant scales 'exit the horizon' during inflation. There are many flat directions which can potentially 
undergo symmetry breaking during inflation, so it is not unlikely that several phase transitions occurred in the 8 
e-folds which is sampled by observations [1] . The observation that the curvature perturbation appears to cut off above 
the scale of the present Hubble radius suggests that (this last period of) inflation may not have lasted much longer 
than the minimum necessary to generate an universe as big as the present Hubble volume. Whereas this raises a 
'naturalness' issue, it is consistent given this state of affairs to consider the possibility that thermal phase transitions 
in flat direction fields ~ 10—15 e-folds after the beginning of inflation leave their mark in the observed scalar density 
perturbation on the microwave sky and in the large-scale distribution of galaxies. 



A. 'Step' model 



The potential for the inflaton (j) 
previously [13]: 



and flat direction fleld (with mass m and /i respectively) is similar to that given 



t < ti, 
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(1) 



Here ti is the time at which the phase transition starts (at i < ti, ^ is trapped at the origin by thermal effects), A is 
the coupling between the and ?/; fields, and 7 is the co-efficient of the leading non-renormalisable operator of order 
n which lifts the potential of the fiat direction field tl^ (and is determined by the nature of the new physics beyond 
the effective field theory description). Note that the quartic coupling above is generated by a term ncfxf)'' ip"^ / Mp in 
the Kahler potential with k ^ 1 (allowed by symmetry near (f> ^ 0), so X = nH^/Mp [H. We have not considered 
non-renormalisable operators cx (/)"/Mp~^ [l^ since we are concerned here with the first ^ 10 — 20 e-folds of inflation 
when (/) is still close to the origin so such operators are then unimportant for its evolution. We are also not addressing 
here the usual '77-problem' in supergravity models, namely that m ^ H due to supersymmetry breaking by the vacuum 
energy driving inflation. We assume that some mechanism suppresses this mass to enable sufficient inflation to occur 
[47I l48j . However this will not be the case in general for the flat direction fields, so one would naturally expect ij ^ H. 
Then the change in the inflaton mass-squared after the phase transition is 



5mi = AS' 



E = 



MS 



A</>^) 



l/(n-2) 



where S is the vev of the global minimum to which ip evolves during inflation. The equations of motion are 



n— 4 



(2) 

(3) 
(4) 



To characterise the comoving curvature perturbation TZ [49j . we employ the gauge-invariant quantity u — —zTZ, 
where z — acp/H, a is the cosmological scale-factor and H is the Hubble parameter during inflation. The Fourier 
components of u satisfy the Klein-Gordon equation of motion: u'^ + (fc^ — z" / z)uk — 0, where the primes indicate 
derivatives with respect to conformal time rj ~ Jdt/a = —1/aH (the last equality holds in de Sitter space). For 

convenience, we use the variable = V2kuk, for which this equation reads: 



+ wi + 



P exp {-2t) - 2-m'^ + Xii; 



2 -2 2Xipip'(j) 



Wk = 0. 



(5) 



where we have used the dimensionless variables: t = Ht, (p = cp/Mp, ip — ip/Mp, m — m/H, X = XMp/H^ and 
k ~ k/Ko, where Kq = aoH and a = aoexp(i). We also define H = H/Mp and Jl = i-i/H for convenience. Note that 
now (and henceforth) the dashes refer to derivatives with respect to t. 
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We start the integration at an initial value of the scale- factor oq, which is taken to be several e- folds of inflation 
before the phase transition occurs. Similarly we choose an initial value (j^o for the inflaton field corresponding to 
several e- folds of inflation before the phase transition occurs (the precise value does not affect our results) . 

We use two parameters to characterise the primordial perturbation spectrum. The first is the amplitude of the 
spectrum on large scales in the slow-roll approximation: 

P,<'>.(JC)'.J^^ ,0) 

\2tt<PoJ 47r2m40g 

The second is ki = cKq where c is a constant chosen so that fci is the position of the step in the spectrum. This can 
be done because a mode with wavenumber k 'exits the horizon' when the coefhcient of Wk in eq.® is zero, hence the 
wavenumber k[ of the mode that exits at the start of the phase transition is 

k[= Kak'i ^Ko\/2 + m^ e*^, (7) 

where ti is just the number of e-folds of inflation after which the phase transition occurs. Since we keep and ti 
fixed, we have k'l oc Kq. For the phase transition to significantly influence the primordial perturbation spectrum after 
the mode with k = k'l exits the horizon requires several further e-folds of inflation, depending on how fast the flat 
direction evolves. The exponential growth of is governed by the value of pt, which we also keep fixed. Therefore the 
position of the step in the spectrum is proportional to kg. We repeat our analysis for integer values of n, the order 
of the non-renormalisable term in the effective field theory description, in the range 12 to 17. We set A = 7 = 1, 
= 0.005, (j)o = 0.01 and Jl^ — 3. The fractional change in the inflaton mass-squared due to the phase transition is 

Am' = — ^ — [ ^ . (8) 

Thus fixing the value of n does not entirely fix Am^, because from eq.(l6]) varying also alters H. A typical 
value is p''^ '-^ 10~^ so the corresponding Hubble parameter is H ^ 3 x 10~^Afp i.e. an inflationary energy scale 
of 2 X lO^** GeV. This is comfortably within the upper hmit of ^ 2 x 10^^ GeV set by the WMAP bound on 
inflationary gravitational waves [50| . Indeed the gravitational wave background is completely negligible for the 'new' 
inflation potential we consider, with the tensor to scalar ratio expected to be: r ~ m^^^ ~ 10~^. 

We varied four parameters describing the homogeneous background cosmology which is taken to be a spatially flat 
Friedman- Robertson- Walker (FRW) universe: the physical baryon density o^b = Slb^i^, the physical cold dark matter 
density ujc = flch^, the ratio 9 of the sound horizon to the angular diameter distance (multiplied by 100), and the 
optical depth r (due to reionisation) to the last scattering surface. The dark energy density is given by f^A = 1 — ^m, 
where fim = + f^b- 

B. 'Bump' model 

Next we consider a multiple inflation model with two successive phase transitions caused by 2 flat direction fields, 
ipi and 'ijj2- The potential is now: 



Vo - \m'^'' - l^il^l + IXicf'i^l + ^i^V"^ (9) 



where ti and t2 are the times at which the first and second phase transitions occur. After i2, for example, the 
equations of motion are 

<^ + 3iJ</)= = (m2 - Ai^A? + AaV^I) (10) 

^2 + 3Hi,2 = -fj = + A20^ - -^ir^-') V'2, (12) 
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- Ai?/;! + >^2i>2 ; \ ■■ • (13) 



and 



The inflaton mass thus changes due to the phase transitions as 

^ rr? - \iY.\ ^ n? - + \2T?2. (14) 

By choosing A2 to be of opposite sign to Ai (possible in the supergravity model), a bump is thus generated in Vtz- 
In the slow-roll approximation the amplitude of the primordial perturbation spectrum first increases from Vn^^^ to 
Vn'^'^ then falls to 'Pn''^\ moving from low to high wavenumbers, where 

■p^m ^ 7^7.(1) = Vn^^^^ (15) 

4^2^Vo' (l~Am2)'' (1- Am^ + Ami)' ■ 

We calculate the curvature perturbation spectrum for this model in a way similar to that for the model with one 
phase transition. As before we specify the homogeneous background cosmology using Wb, and r and, anticipating 
the discussion later, the fraction of dark matter in the form of neutrinos fi, = 57^ /fid where the total dark matter 
density is U,^ = Q.^, + The primordial perturbation spectrum is parameterised in a similar way to that of the first 
model using V-jz-^^ and ki = cKq, which is now the approximate position of the first step in the spectrum. A measure 
of the position of the second step is given by the third parameter 

fca = fcie(*'^-*i\ (16) 

where the exponent is just the number of Hubble times after the beginning of the first phase transition when the 
second one starts. Initially we examined models with ni — 12 and 712 = 13, and then we let V-r}-'^'^ and Vn-"^^ vary 
freely. We set Ai, A2, 71 and 72 all equal to unity, rh? = 0.005, (pQ = 0.01 and 'jl\ = 'jl\ = ?, throughout. 



III. THE DATA SETS 

We fit to the WMAP 3-year Q temperature-temperature (TT), temperature-electric polarisation (TE), and electric- 
electric polarisation (EE) spectra alone as we wish to avoid possible systematic problems associated with combining 
other CMB data sets. We also fit the linear matter power spectrum Pni{k) to the SDSS measurement of the real space 
galaxy power spectrum Vg{k) [33| . The two spectra are taken to be related through a scale-independent bias factor 
^SDSS so that Vm{k) = ^sDSs'^g(^)- This bias is expected to be close to unity and we analytically marginalise over it 
using a flat prior. 

We also fit our models to the redshift two-point correlation function of the SDSS LRG sample, which was obtained 
assuming a fiducial cosmological model with ~ 0.3, flA — 0.7 and h = 0.7 The correlation function is 

dependent on the choice of background cosmology so we rescale it appropriately as explained below, in order to 
confront it with the other cosmological models we consider. Again we consider a linear bias 6lrG: and determine it 
from the fit itself. 



A. WMAP 



To date the most accurate observations of the CMB angular power spectra, both TT and TE, have been made by 
the WMAP satellite using 20 differential radiometers arranged in 10 'differencing assemblies' at 5 different frequency 
bands between 23 and 94 GHz. Each differencing assembly produces a statistically independent stream of time-ordered 
data. Full sky maps were generated from the calibrated data using iterative algorithms. The use of different frequency 
channels enable astrophysical foreground signals to be removed and the angular power spectra were calculated using 
both 'pseudo-Cf' and maximum likelihood based methods [6|. We use the 3-year data release of the TT, TE and EE 
spectra and also the code for calculating the likelihood (incorporating the covariance matrix), made publicly available 
by the WMAP team.^ 



^ http://lambda.gsfc.nasa.gov/ 
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Parameter 


Lower limit 


Upper limit 




0.005 


0.1 




0.01 


0.99 


e 


0.5 


10.0 


T 


0.01 


0.8 


In (lOlOp7^(°>) 


0.01 


6.0 


10*fci/Mpc'i 


0.01 


600.0 


6lrg 


1.0 


4.0 


h 


0.4 


1.0 


Age/Gyr 


10.0 


20.0 



TABLE I: The priors adopted on the input parameters of the step model, as well as on the derived parameters: the Hubble 
constant and the age of the Universe. 



Parameter 


Lower limit 


Upper limit 




0.005 


0.1 


e 


0.5 


10.0 


T 


0.01 


0.8 




0.01 


0.3 


In (lO^opT^f")) 


0.01 


6.0 


6lrg 


1.0 


4.0 


h 


0.1 


1.0 


Age/Gyr 


10.0 


20.0 



TABLE IL The priors adopted on the input parameters of the bump model with n\ — 12, ni = 13, as well as on the derived 
parameters: the Hubble constant and the age of the Universe. The priors set on the position of the steps are given in Table HITl 

Our special interest is in the fact that apart from the quadrupole there are several other multipoles for which the 
binned power hes outside the la cosmic variance error. These are associated with small, sharp features in the TT 
spectrum termed glitches, e.g. at ^ = 22, 40 and 120 in both the 1-year and 3-year data release. The pseudo-Q 
method produces correlated estimates for neighbouring C^'s which means it is difficult to judge the goodncss-of-fit of 
a model by eye. The contribution to the P^r multipole for the best fit ACDM model is shown in Fig. 17 of ref.Q, 
where it is apparent that the excess originates largely from the multipoles at ^ < 120. 

B. SDSS 

The Sloan Digital Sky Survey (SDSS) consists of a 5 band imaging survey and a spectroscopic survey covering 
a large fraction of the sky, carried out together using a 2.5 m wide- field telescope. For the spectroscopic survey 
targets are chosen from the imaging survey to produce three catalogues: the main galaxy sample, the luminous, red 
galaxy (LRG) sample and the quasar sample, which extend out to redshifts ol z 0.3, 0.5 and 5 respectively. The 
matter power spectrum (up to a multiplicative bias factor &sdss) has been measured using about 2 x 10^ galaxies 
which form the majority of the main galaxy sample from the second SDSS Data Release [33|. After correcting for 





WMAP 


+SDSS 


+LRG 


+SDSS+LRG 


10''fci/Mpc~^ 


600.0 


600.0 


600.0 


600.0 


lO^fe/Mpc"^ 


1100.0 


1100.0 


600.0 


600.0 



TABLE HL The prior upper limits set for the position of the steps in the spectrum for the 'bump' model. In each case the 
lower limit is 10"'' Mpc"^ 
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redshift-space distortion due to galaxy peculiar velocities, the galaxy density field was expanded in terms of Karhunen- 
Loeve eigenmodes (which is more convenient than a Fourier expansion for surveys with complex geometries) and the 
power spectrum was estimated using a quadratic estimator. We fit to the 19 k-band measurements in the range 
0.01 < k < 0.2 Mpc^^, as the fluctuations turn non-linear on smaller scales. 



C. LRG 



The large effective volume of the SDSS LRG sample allowed an unambiguous detection of the BAO peak in the 
two-point correlation function of galaxy clustering [3a]- The redshifts of 5 x lO'' LRG's were translated into comoving 
coordinates assuming a fiducial Qjy — 0.7, fl^i = 0.3 cosmology. The correlation function was then found using the 
Land-Szalay estimator. This can be rescaled for other cosmological models as follows. 

An object at redshift z with an angular size of A0 and extending over a redshift interval Az has dimensions in 
redshift space xy and x±, parallel and perpendicular to the line-of-sight, given by 



^11 



Az 



a;_L = {1 + z)Da{z)A9. 



(17) 



Thus the volume in redshift space is cx D\{z)/H{z), where I?a is the angular diameter distance. The distances 
between galaxies in redshift space, hence the scale of features in the redshift correlation function, are proportional to 
the cube root of the volume occupied by the galaxies in redshift space. Therefore the multiplicative rescaling factor 
between the scale of the BAO peak in the fiducial model and that in another model is just: 



Dl{z)HUz) 



DlU^)H{z) 



(18) 



i.e. Xpcak = 7rcs X 2:pcak,fid- After rcscaliug in this way for the cosmological models we consider, the 
the fit of ^ to the data points was found using the appropriate covariance matrix. 



statistic for 



IV. CALCULATION OF THE OBSERVABLES 

Since MCMC analysis of inflation requires the evaluation of the observables for typically tens of thousands of 
different choices of cosmological parameters, it is important that each calculation be as fast as possible. 

The inflaton mass is constant before and after the phase transition(s), so the equations of motion for cj) and have 
the solutions 



= e-*/2 



cie*'-+ H-C2e*'--, (19) 
C3H(,i) (fee-*') + oRf^ (^e-*)l . (20) 



Here ci to C4 are constants, hI.^'' and hI?-* are Hankel functions of order v — y j — and r± — — 3/2± z^. We evolve 
Tp and Wk numerically through the phase transition(s), matching to eas. (|19p and (j20p . The slow- roll approximation 
is applicable well before the phase transition(s), so the initial values are related as 0q = {m? /'i)(t>Q. We begin the 
integration of the Klein-Gordon equation at igtart when the condition 



Id^ 

z drf 



= (21) 



is satisfied. The arbitrary constant e is chosen to be 5 x 10^^ which is sufficiently small for V-ji to be independent of 
the precise value of tstart, yet lar ge e nough for the numerical integration to be feasible. For the Bunch-Davies vacuum, 
the initial conditions for Wk are [49| 

U'fe(?'start) = 1, w'k(tstait) = -ifccxp {-tstart) ■ (22) 



^ http://cmb.as.arizona.edu/ eisenste/ 
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The amplitude of the primordial perturbation spectrum on large scales was given in eq.®; on smaller scales it is 

evaluated when the mode has crossed well outside the horizon. Rather than calculating V-r, in this way for every 
wavenumber, we use cubic spline interpolation between 0(100) wavenumbers ki to save time. These sample the 
spectrum more densely where its curvature is higher, in order to minimise the interpolation error. The ki values are 
found each time using the adaptive sampling algorithm described in Appendix |^ 

Taking the interpolated primordial perturbation spectrum as input, we use a modified version of the cosmological 
Boltzmann code CAMB [51| to evaluate the CMB TT, TE and EE spectra, as well as the matter power spectrum.^ 

We find the two-point correlation function using the following procedure: for 10~^ < k < 2/iMpc~^ the linear 
matter power spectrum Pm{k) is obtained using CAMB, while outside this range baryon oscillations are negligible, 
and so to calculate the matter power spectrum Pm(k) for 10"^ < k < 10^^ hMpc^^ and 2 < k < 10^ /iMpc~^ we use 
the no-baryon transfer function fitting formula of ref. [53] normalised to the CAMB transfer function. This significantly 
increases the speed of the computation. The 'Halofit' procedure [53| is then used to apply the corrections from non- 
linear evolution to produce the non-linear matter power spectrum (k). The LRG power spectrum is given by 
PhRoik) = ^LRG^m'"(^)- "^^^ ^^^^ spacc Correlation function ^c{x) is then calculated by taking the Fourier transform. 
The integral is performed over a finite range of k using the FFTLog code [sj which takes the fast Fourier transform of 
a series of logarithmically spaced points. The wavenumber range must be broad enough to prevent ringing and aliasing 
effects from the EFT over the interval of x for which we require ^c- We find the range 10~^ < k < 10^ ZiMpc"^ to be 
sufficient. Finally, the angle averaged redshift correlation function is found from the real space correlation function 
using 

a^)-{i + lf+lf)u^), /^^' (24) 

which corrects for redshift space distortion effects [ssl . [s^ . 



V. RESULTS 



We calculate the mean values of the marginalised cosmological parameters together with their 68% confidence 
limits, for the four combinations of data sets: WMAP, WMAP + SDSS, WMAP + LRG, WMAP + SDSS + LRG. 
In Table HVl we show this for the ACDM model with a Harrison- Zeldovich spectrum. Adding an overall 'tilt' to the 
spectrum creates a noticeable improvement (by Ax^ ^ 7 — 9) as is seen from Table |V] which gives results for the 
ACDM model with a power-law scalar spectral index of Ug ~ 0.95 (pivot point k = 0.05 Mpc^^). The goodness-of- 
fit statistic is calculated using the WMAP likelihood code, supplemented by our estimate for the LRG fits obtained 
using the covariance matrix used in ref. 38]. 

In order to quantify the desirable compromise between improving the fit and adding new parameters, it has become 
customary to use the "Akaike information criterion" defined as AIC = — 21n£,nax + 2iV [s^], where £max is the 
maximum likelihood and N the number of parameters. Although not a substitute for a full Bayesian evidence 
calculation (which is beyond the scope of this work), this is a commonly used guide for judging whether additional 
parameters are justified — the model with the minimum AIC value is in some sense preferred (58| . According to 
this criterion, the power-law ACDM model is preferred over the scale-invariant ACDM model, hence we consider the 
former to be the benchmark against which our models should be judged. However we wish to emphasise that we have 
physical motivation for the non-standard primordial spectra we consider. We are not performing a pure parameter 
fitting exercise for which criteria like the AIC might be more relevant. 



A. ACDM 'step' model 

In Tables IVII to IIXI we show that adding a step to the spectrum improves the fit by Ax^ 1 — 2 relative to the 
scale-invariant model, but not as much as addition of a tilt. Since the model has 2 additional parameters — the 
position of the step ki and n (which determines Am^) — Aaic is positive relative to the power-law ACDM model. 
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TABLE IV: lu constraints on the marginalised cosmological parameters for a scale-invariant ACDM model. The 6 parameters 
in the upper part of the Table are varied by CosmoMC, while those in the lower part are derived quantities. The of the fit 
is given, as is the Akaike information criterion relative to the power-law ACDM model in Table FVl 
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TABLE V: la constraints on the marginalised cosmological parameters for the power-law ACDM model. The 6 parameters in 
the upper part of the Table are varied by CosmoMC, while those in the lower part are derived quantities. The of the fit is 
given as is the relative Akaike information criterion (normalised to be zero for this model). 

We have not shown the results for the ACDM step model allowing an overall tilt in the spectrum; although the 
fits do improve again (by Ax^ ~ 1 — 2) relative to the power-law ACDM model, they are not favoured by the AIC 
because of the 2 additional parameters ki and n. The values of the background cosmological parameters do not change 
significantly with n either, and remain similar to those of the power-law ACDM model. However, the values of both 
Vn''^^ and fci fall with increasing n, as seen in the 1-D likelihood distributions for the parameters shown in FiglTJ 

To understand this, recall that the amplitude of the primordial perturbation spectrum is well constrained by the 
TT spectrum on medium scales, but is uncertain on large scales due to cosmic variance. Therefore the data can 
accommodate a large step in V-r provided that the top of the step is at the right height to match the TT spectrum. 
Consequently V-r}-^^ must vary with the infiaton mass change in such a way as to leave the primordial perturbation 
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TABLE VI: la constraints on the marginalised cosmological parameters for the ACDM step model using WMAP data alone. 
The 6 parameters in the upper part of the Table are varied by CosmoMC, while those in the lower part are derived quantities. 
The of the fit is given, as is the Akaike information criterion relative to the power-law ACDM model in Table FVl (taking into 
account that n is also a parameter). 
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TABLE VII: la constraints on the marginalised cosmological parameters for the ACDM step model using WMAP + SDSS 
data. The 6 parameters in the upper part of the Table are varied by CosmoMC, while those in the lower part are derived 
quantities. The of the fit is given, as is the Akaike information criterion relative to the power- law ACDM model in Table FVl 
(taking into account that n is also a parameter). 

spectrum invariant on smaller scales. This is seen in Figl^] which shows the primordial spectra for the best-fit models. 

The position of the step in Vqi^k) for n = 14 means that the associated TT and TE spectra are lower on large 
scales than those of the best-fit power-law model — see Figsl3]and|31 For n = 15 to n = 17 the lower plateau of the 
step is too far outside the Hubble radius to suppress the TT and TE spectra. However the primordial spectrum has 
a prominent first peak at fc ~ 6 x 10^** h Mpc~^. This introduces a maximum in the TT spectrum centred on £ = 4 
which fits the excess power seen there by WMAP. There is a corresponding peak at low £ in the TE spectrum. The 
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TABLE VIII: Icr constraints on the marginalised cosmological parameters for the ACDM step model using WMAP + LRG 
data. The 7 parameters in the upper part of the Table are varied by CosmoMC, while those in the lower part are derived 
quantities. The of the fit is given, as is the Akaike information criterion relative to the power- law ACDM model in Table FVl 
(taking into account that n is also a parameter). 

peak at fc ~ 1.5 x lO^'^ /iMpc^' in the primordial spectrum (for n = 14 to n = 17) increases the amplitude of the XT 
spectrum around ^ = 15. However the glitches in the TT spectrum at = 22 and £ — 40 are too sharp to match the 
oscillations produced by the mechanism considered here. 

As seen in FiglBl the SDSS galaxy power spectrum is well matched in all cases but the features in the spectrum due 
to the phase transition appear far outside the scales probed by redshift surveys (see FiglT]). Moreover a step in the 
matter power spectrum does not significantly alter the two-point correlation function in a ACDM universe as shown 
in FigO Just as for the concordance power-law ACDM model however, the amplitude of the predicted BAO peak is 
too low by a factor of ^ 2, although its predicted position does match the data (37| . 

We conclude therefore that the observations favour a 'tilted' primordial spectrum over a scale-invariant one, as 
has been emphasised already by the WMAP team However the spectrum is not required to be scale-free. Very 
different primordial spectra (as shown in Figl2|) provide just as good a fit to the data, although this is admittedly 
penalised by the Akaike information criterion taking into account the 2 additional parameters characterising the step 
in the spectrum. This motivates us to ask whether even the best-fit cosmological model might be altered if a more 
radical departure from a scale- free spectrum is considered. For example in multiple inflation [Sj , two phase transitions 
can occur in rapid succession resulting in a 'bump' in the primordial spectrum. Such a feature has been advocated 
earlier on empirical grounds for fitting a ACDM model [5§|. 

B. CHDM 'bump' model 

We consider now a primordial spectrum with a bump at fc ~ 2 x 10~'^/iMpc~' generated by 2 successive phase 
transitions which cause an upward step followed by a slightly larger downward step in the amplitude of the primordial 
perturbation, as shown in FiglS] This boosts the amplitude of the TT spectrum on the left of the first acoustic peak 
but suppresses the second and third peaks. This is all that is necessary to fit the WMAP data to an Einstein-de 
Sitter cosmology as seen in FiglTU] Since there is no late ISW effect, the amplitude is smaller on large scales than in 
an universe with dark energy. The fits to the TE and EE spectra are also good as shown in FigUT] and [T^l 

By adding a hot dark matter component in the form of massive neutrinos, the amplitude of the matter power 
spectrum on small scales is suppressed relative to a pure CDM model and a good fit obtained to SDSS data (see 
FigfT3|. As noted earlier [H], this suppresses erg so as to provide better agreement with the value deduced from 
clusters and weak lensing. A further suppression occurs due to the downward step in the primordial spectrum. 
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TABLE IX: Icr constraints on the marginalised cosmological parameters for the ACDM step model using WMAP + SDSS + 
LRG data. The 7 parameters in the upper part of the Table are varied by CosmoMC, while those in the lower part are derived 
quantities. The of the fit is given, as is the Akaike information criterion relative to the power- law ACDM model in Table FVl 
(taking into account that n is also a parameter). 

The constraints on the marginalised cosmological parameters are given in Table |X] with the order of the non- 
renormalisable terms set to ni — 12 and n2 — 13 (which we nevertheless count as 2 additional parameters). We also 
allow ni and n2 to vary continuously (but find no further improvement), with the results shown in Table IXll and as 1-D 
likelihood distributions in Fig[T4l The fit of this model to the WMAP data and the SDSS matter power spectrum is 
just as good as that of the pow er-law ACDM model, as indicated by the values and the vanishing Aaic- However, 
as has been noted already [S^, the CHDM model does not fit the LRG two-point correlation function well. This is 
because the BAO peak corresponds to the comoving sound horizon at baryon decoupling; the latter is a function of 
rim and is too low in an Einstein-de Sitter universe as seen in Fig llSI 
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FIG. 1: 1-dimensional likelihood distributions of the marginalised cosmological parameters for the ACDM step model using the 
WMAP + SDSS + LRG data. 



VI. PARAMETER DEGENERACIES 



If two different parameters have similar effects on an observable then they will be negatively correlated in a 2-D 
likelihood plot since the effect of increasing the value of one of them can be undone by decreasing the value of the 
other. On the other hand, if two parameters have opposite effects then increasing the value of one can compensate 
for increasing the value of the other and the parameters will be positively correlated. Such correlations are known as 
parameter degeneracies and they limit the constraints which the data can place on the parameter values. 



A. ACDM 'step' model degeneracies 



The results indicate a strong parameter degeneracy amongst the models with one phase transition. Moving in 
parameter space along the direction of the degeneracy, Am^ increases, while Vn-^^ and ki both decrease. This 
arises because the amplitude of the TT spectrum on the scale of the acoustic peaks is governed by the combination 
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FIG. 2: The best-fit primordial perturbation spectra for the ACDM step modeL Note the suppression of power at the 
wavenumber coresponding to the present Hubble radius: l/o — 3 x 10~* /iMpc"^. 
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FIG. 3: The best-fit TT spectra for the ACDM step model, with WMAP data. 
■PTj^^-'e"^'^, where V-rS^^ is the amplitude of the primordial spectrum after the phase transition: 
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The fit to WMAP data thus requires V-ji-^^ oc e^"^. Consequently r at fixed Am? is positively correlated with Pt?.'-*'' 
(reflecting the well-known degeneracy between the optical depth and the normalisation of the TT spectrum). However 

r falls going from n = 15 to n = 17 because decreasing r reduces P^\ while increasing Am^ raises P^'^ ■ 

The value oiVu'^^ is constrained after marginalising over r, which means that the relationship between Am^ and 
is fitted weU by 



p4") = (1 - Am^) 



(26) 
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FIG. 4: The best-fit TE spectra for the ACDM step model, with WMAP data. 
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FIG. 5: The best-fit EE spectra for the ACDM step model, with WMAP data. 



where P™ is a constant (being the amplitude of a scale- invariant Harrison-Zeldovich primordial spectrum) . This is 
illustrated in FiglTH] using results from Table IIXI 

The parameter fci is negatively correlated with Am^ because the higher cosmic variance at low I allows the data to 
accommodate more prominent features in the primordial spectrum at smaller k. The likelihood distribution of ki is 
strongly non-Gaussian, as shown in FigU] Each distribution has maxima wherever the oscillations in the corresponding 
TT spectrum due to the phase transition match the glitches in the WMAP measurements. 



B. CHDM 'bump' model degeneracies 



The CHDM bump model has several parameter degeneracies which are illustrated in FiglTT] There is a positive 
correlation between the optical depth and the amplitude of the primordial spectrum on medium and small scales. 



16 



10* 



10^ 




ACDM step, n = 14 

ACDM step, n = 15 

ACDM step, n = 16 

ACDM step. 71 = 17 

ACDM power-law 



10- 



10- 



k (ft Mpc-'J 

FIG. 6: The best-fit matter power spectra for the ACDM step model, with SDSS data. 
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FIG. 7: Oscillations in the matter power spectra on very large scales in the ACDM step model. 



The correlation between r and is weaker because reionisation does not damp the TT spectrum on large scales. 

Increasing the baryon density increases the height of the first acoustic peak relative to the second one [6(1] , while 
increasing Am| increases the size of the second step in the primordial spectrum and so boosts the height of the first 
peak relative to the other peaks. Thus ujh and Ato^ ^re negatively correlated since both parameters have the same 
effect on the height of the first peak relative to the second [60| ■ 

The bump in the primordial spectrum becomes broader when ki is increased, which increases the height of the 
second acoustic peak relative to the third. Reducing the density of CDM has a similar effect on the peak heights, 
hence ki and ujc are positively correlated. The baryon and CDM densities are also positively correlated because 
increasing Wb increases the height of the first acoustic peak, while increasing lUc has the opposite effect (60| . 

The error bars of the WMAP data are smallest in the region of the first acoustic peak, so P^"^ is more tightly 
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FIG. 8: The best-fit two-point galaxy correlation functions in the ACDM step model, with LRG data; the spatial scales have 
been shifted using eq.((T8} to match the same ACDM cosmology and enable comparison. 
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FIG. 9: The primordial perturbation spectrum for the CHDM bump model with ni — 12 and n2 — 13, compared to the ACDM 
power-law model with Us — 0.95. 



constrained than Pj^ or Pj^ , as is apparent from Table IXII Consequently the parameters P-j^ and Ami satisfy 

P4°)=PH^(1-A™?)^ (27) 

which is similar to eq. (|26|) . 

From the last equality of eq. p3|) it is apparent that increasing Am^ reduces Vn^'^K hence these parameters are 
anti-correlated, as seen in the final panel of FigfTTl 
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FIG. 10: The best-fit TT spectrum for tlie CHDM bump model, with WMAP data. 
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FIG. 11: The best fit TE spectrum for the CHDM bump model, with WMAP data. 



VII. CONCLUSIONS 



Anisotropics in the CMB and correlations in the large-scale distribution of galaxies reflect the primordial perturba- 
tions, presumably from inflation, but convoluted with the a priori unknown effects of their evolution in matter. Thus 
the values of the parameters of the assumed cosmological model are necessarily uncertain due to our limited knowledge 
of both the dark matter content of the universe, as well as of the nature of the primordial perturbations. Nevertheless, 
by combining various lines of evidence, it has been possible to make a strong case for the concordance ACDM model, 
based on a FRW cosmology [6l|. This model appears deceptively simple, but invokes vacuum energy at an unnatural 
scale, Pa ~ (lO^'^'^Afp)"', as the dominant constituent of the universe. There is no satisfactory understanding of this 
from fundamental theoretical considerations j6^ . [63j . It is therefore worth reexamining whether the precision WMAP 
data, which has played a key role in the general acceptance of the ACDM model, cannot be otherwise interpreted. 

In this work we have focussed on the fact that our present understanding of the physics of inflation is rather primitive 
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FIG. 12: The best fit EE spectrum for the CHDM bump model, with WMAP data. 
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FIG. 13: The best-fit matter power spectra for the CHDM bump model, with SDSS data. 



and there is no compelling reason why the generated primordial adiabatic scalar density perturbation should be close 
to the scale-invariant Harrison-Zeldovich form, as is usually assumed. This is indeed what is expected in the simplest 
toy model of an inflaton field evolving slowly down a suitably flat potential, but attempts to realise this in a physical 
theory such as supergravity or string/M-theory are plagued with difficulties. Such considerations strongly suggest 
moreover that even if just one scalar field comes to dominate the energy density and drives inflation, the non-trivial 
dynamics of other scalar fields in the inflating universe can affect its slow-roll evolution and thus create features in the 
primordial spectrum. Indeed there is some indication for such spectral features in the WMAP observations. Whether 
they are real rather than systematic effects will be tested further in the WMAP 5-yr analysis and by the forthcoming 
Planck mission [64{ . It has also been noted that significant non-Gaussianity would be generated when there are sharp 
features in the inflaton potential [65]. This can be computed in our model and would constitute an independent test. 

We have demonstrated that such spectral features can have a major impact on cosmological parameter estimation, 
in that the WMAP data can be fitted without requiring any dark energy, if there happens to be a mild bump 
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TABLE X: la constraints on the marginalised cosmological parameters for the CHDM bump model, with the changes in the 
inflaton mass set by choosing ni = 12, 712 = 13 (which are still counted as 2 additional free parameters). The 8 parameters in 
the upper part of the Table are varied by CosmoMC, while those in the lower part are derived quantities. The of the fit is 
given, as is the Akaike information criterion relative to the power- law ACDM model in Table [V] 

in the primordial spectrum on spatial scales of order the Hubble radius at the last scattering surface. This has 
been modelled here using the supergravity-based multiple inflation model Q with reasonable choices for the model 
parameters. Although other possibilities for generating spectral features have been investigated, our framework has 
the advantage of being based on effective field theory and thus calculable. 

We have also departed from earlier work on parameter estimation in considering possible values of the Hubble 
constant well below the HKP value ^33| which is usually input as a prior. Moreover the CHDM 'bump' model cannot 
fit the luminosity distance measurements of SN la, since the supernovae are fainter than predicted in a homogeneous 
Einstein-de Sitter cosmology. However, it has been noted that deep measurements using physical methods yield a 
lower value of h than that measured locally [13, as would be the case if we are located in an underdense void which 
is expanding faster than the average [6^. There are ongoing attempts to test whether we are indeed inside such a 
'Hubble bubble' [67l.[68|. This is enormously important since such an inhomogeneous LTB cosmology can in principle 
address the failures of our FRW framework, with regard to the SN la Hubble diagram and the B AO peak [33] ■ 

Moreover we have shown that the observed power spectrum of galaxy clustering can also be accounted for by 
modifying the nature of the dark matter, in particular allowing for a significant component due to massive neutrinos. 
Whether neutrinos do have the required mass of ^ 0.5 eV will be definitively tested in the forthcoming KATRIN 
experiment (69j . The position of the baryon peak seen in the galaxy correlation function may also be matched in the 
LTB framework, with a low global Hubble constant [39*1 . 

Regardless of whether this model is the right description of the physical universe, we wish to emphasise that it 
provides a good fit to the WMAP data. Thus for progress to be made in pinning down the cosmological parameters 
and extending our understanding of inflation using precision CMB data, it is clearly necessary for a broader analysis 
framework to be adopted than has been the practice so far. 
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TABLE XI: la constraints on the marginalised cosmological parameters for the CHDM bump model, with the changes in the 
inflaton mass allowed to vary freely. The 10 parameters in the upper part of the Table are varied by CosmoMC, while those in 
the lower part are derived quantities. The oi the fit is given, as is the Akaike information criterion relative to the power-law 
ACDM model in Table |Vl 
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APPENDIX A: ADAPTIVE SAMPLING ALGORITHM 

For plotting or interpolating a function y — f{x) over an interval (a^min, a^max) it is useful to have an algorithm that 
generates a set of points {xi, f (xi)}, where x^+i > Xi, such that the density of the samples Xi in the abscissa increases 
with d^f/dx^. One such algorithm is listed here. It works by ensuring that the value of the function at the mid-point 
of each interval {xi,Xi+2) is within some tolerance e of the value there linearly interpolated from the end-points of 
the interval. 

1. Form two lists {xi, X2t ■ ■ ,xn} and {j/i, 2/2, ■ • ■ , Un} where the x'^s are equally spaced, 

— Xnun 1 ^ -^maxi ^ — 1:2,..., A'^, (.^1) 

and yi = f (x,,). 

2. Start with the first interval {xi,X2). 

3. For the interval {xi,Xi+i) under consideration calculate Xm = {xi + Xi+i)/2 and ym = f{xm)- Insert Xm 
between Xi and Xi+i in the list of x values and ym between yi and yi+i in the list of y values. 

4. If \ym — (yi + 2/i+i) /2| > e repeat step 3 with the new interval (xi, Xm)- 
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FIG. 14: 1-dimensional likelihood distributions of the marginalised cosmological parameters for the CHDM bump model using 
the WMAP, SDSS and LRG data. 

5. Otherwise repeat step 3 with the next interval (xi+i, unless x^+i — Xmax, in which case finish. 

This algorithm requires a minimum number of function evaluations and none are wasted. In our work we set a; = In fc 
and achieve satisfactory results with = 20 and e = 10^^^, as illustrated in FigfTSl 

APPENDIX B: MCMC LIKELIHOOD ANALYSIS 



We use the Monte Carlo Markov Chain (MCMC) approach to cosmological parameter estimation. It is a method 
for drawing samples from the posterior distribution P (roldata) of the parameters tu, given the data. According to 
Bayes' theorem, the posterior distribution is given by: 



F(Tj7|data) = 



C {data\zu) P (v^) 
J C (dataj-ro-) P (-cc?) dzu ' 



(Bl) 



where C (dataj-ccr) is the likelihood of the parameters and P {txt) is their prior distribution. We use flat priors on all 
the parameters listed in Table HI With the advent of massively parallel supercomputers which can run several chains 
simultaneously MCMC has become the standard tool for model fitting. At best, the number of likelihood evaluations 
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FIG. 15: The best-fit two-point galaxy correlation functions in the CHDM bump model, with LRG data. The best-fit ACDM 
power-law model, appropriately rescaled, is shown for comparison. 
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FIG. 16: The degeneracy between the large-scale amplitude of the primordial spectrum and the fractional change in the inflaton 
mass according to ea.(|25p. with the 'data' points talcen from Table HXl 



required increases linearly with the number of parameters, n. This is much slower than traditional grid based methods 
which need m" evaluations where m is the number of grid points in each parameter. 

We use a modified version of the CosmoMC software package [70].^ CosmoMC contains an implementation of the 
Metropolis-Hastings algorithm for generating Markov chains, which runs as follows: let the vector w'*) be the current 
location of the chain in parameter space. A candidate -cct for the next step in the chain is chosen from some proposal 

distribution q zu'-^^ ) . The proposed location ^ may differ from -ccrf*) by either some or all of the parameter values. 



^ http:/ /cosmologist.info/cosmomc/ 
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FIG. 17: Degeneracies between the marginalised cosmological parameters for the CHDM bump model. 
Then -cct(*+^) — to' with probabihty 

a ro*-*-'^ = mil 



q (t^(^) 


zu^ C ^data 




q (jzr 





(B2) 



and the step is said to be 'accepted', otherwise = ro'*) and the step is 'rejected'. A chain that is built up 

from many such steps will be Markovian as the proposal distribution q ^ro zu^^^^ is a function of -ccr'*) only and not 
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FIG. 18: An example result from our adaptive sampling algorithm which samples the power spectrum more finely where the 
curvature is higher. 



earlier steps. When q (^Tzr 



q (-cl7 — ■cu'^*'), as in the case of CosmoMC, the chain executes a random walk in 



parameter space, apart from the rejected steps. Since a satisfies the reversibility condition, 



t(0 



(B3) 



It can be shown by construction that C {vc) will be the equilibrium distribution of the chain in nearly all practical 
situations. Thus when equilibrium is reached the chain represents a random walk where proportionally more time 
is spent in regions of higher likelihood. This occurs after a transient period known as the 'burn in' where the steps 
are highly correlated with the starting point (and are usually discarded). Once the chain is equilibrated, it must run 
long enough to explore all of the relevent parameter space after which the chain is said to be 'mixed'. A chain which 
moves rapidly through parameter space so that this happens quickly is said to have 'good mixing'. Shot noise in the 
statistics inferred from the chain is also reduced by running longer. 

When the statistics accurately reflect the posterior distribution and the chain can be used reliably for parameter 
estimation, it is said to have 'converged'. In practice, convergence is often taken to have occurred when the results 
derived from a chain are independent of its starting point. Many methods of diagnosing convergence have been 
proposed. CosmoMC uses one recommended in Tei.Jl] which involves running several chains with widely dispersed 
starting points. (For a different technique used in the context of cosmological parameter estimation see ref.[7^.) A 
comparison is made between the within-chain variance of the chains, which is the mean of the variance of each of the 
chains, and the between-chain variance, which is equal to the length of the chains multiplied by the variance of the 
means of the chains (in the case where all chains have the same length) . The idea is that they should be the same 
after convergence, as all the chains exhibit the same behaviour and the dependence on the starting position has been 
lost. Ref.[71] present a quantity R containing the ratio of the within-chain and the between-chain variances such that 
at convergence i? ~ 1. In our work we assume convergence to have occurred when R = 1.02. 

The speed with which a chain converges is strongly dependent on the proposal distribution. If too many large steps 
are proposed, the acceptance rate will be low as most steps will be away from the region of high likelihood, and the 
chain will be stationary for long periods of time. If too many small steps are proposed, most will be accepted but the 
chain will move slowly through parameter space. In either case the chains will mix slowly and the steps will be highly 
correlated. We use the results of preliminary MCMC runs to fine-tune the width of the proposal distribution before 
proceeding with the final run. 

Cosmological parameters in CosmoMC are divided into two classes: 'fast' and 'slow'. Slow parameters are those 
upon which the transfer function is dependent. All others are fast parameters, including those which govern the 
primordial power spectrum or correspond to calibration uncertainties in the data. Steps in which any of the slow 
parameters change take a relatively long time to compute, because the complicated Einstein-Boltzmann equations 
for the transfer function must be solved. Once the transfer function is known however, steps in which only the fast 
parameters change can be taken quickly as CAME uses linear perturbation theory. Although the calculation of the 
primordial power spectrum takes longer than usual in our modified version of CAMB, a fast step is still quicker than 
a slow one, and we retain the fast and slow division of parameters. CosmoMC exploits this split by alternately taking 
fast and slow steps, which allows a more rapid exploration of the parameter space than slow steps alone. 
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In the absence of a covariance matrix for the parameters, CosmoMC chooses at random a basis in the slow parameter 
subspace. It then proposes in turn a step in the direction of each basis vector with every subsequent slow step. Once 
all the basis vectors have been used in this way, CosmoMC chooses a new random set and repeats the cycle. The 
length of the step is determined by the proposal distribution, which in CosmoMC is based upon a two dimensional 
radial Gaussian function mixed with an exponential. Fast steps are taken in a similar way, by cycling through random 
basis vectors in the fast subspace. This is done to reduce the risk of the chain doubling back on itself. 

When a covariance matrix is available, CosmoMC takes slow steps by cycling through random bases in the Uaiow 
dimensional subspace spanned by the Usiow largest eigenvectors of the covariance matrix. In this way fast parameters 
that are correlated with slow ones are also changed during a slow step, in the direction of the degeneracies, which 
increases the mobility of the chain. Fast steps are just made in the fast subspace as before, and together the fast and 
slow steps can traverse the whole parameter space. 

Given N samples -ca^*) drawn from P {■cj \ data) the best estimate for the distribution is formally 



Pn data) = ^ " "^^'0 ' ^^^^ 



However, rather than studying the full function P(-ccT|data) it is often easier to interpret marginalised probability 
distributions obtained by integrating over a subset of the cosmological parameters, 

/n 
P (n7i , 1372 , . . . , n7„ I data) Y\ dwi. (B5) 

z— m+l 

From eq. (|B4p the probability given the data that a parameter lies in a particular interval is proportional to the number 
of samples that fall into the interval. Thus the marginalised probability distributions can be found using histograms 
of the samples. For 1-D marginalised probability distributions confidence intervals are frequently quoted such that a 
fraction 1 — a of the samples fall within the interval, while a/2 lie higher and a/2 lie lower. 
The expectation value of a parameter given the data is 

(ci7a| data) = J ccjaP (tUa I data)d-ci7a. (B6) 

Using eas. (|B4l) and (jB5|) it is approximated by 

1 ^ 

(tnJdata)^ = -^7i7W. (B7) 
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